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ABSTRACT 

In two recent papers the mesoscale model Meso-NH, joint with the Astro-Meso-NH 
package, has been validated at Dome C, Antarctica, for the characterization of the 
optical turbulence. It has been shown that the meteorological parameters (tempera- 
ture and wind speed, from which the optical turbulence depends on) as well as the 
profiles above Dome C were correctly statistically reproduced. The three most impor- 
tant derived parameters that characterize the optical turbulence above the internal 
antarctic plateau: the surface layer thickness, the seeing in the free-atmosphere and in 
the total atmosphere showed to be in a very good agreement with observations. Vali- 
dation of has been performed using all the measurements of the optical turbulence 
vertical distribution obtained in winter so far. In this paper, in order to investigate the 
ability of the model to discriminate between different turbulence conditions for site 
testing, we extend the study to two other potential astronomical sites in Antarctica: 
Dome A and South Pole, which we expect to be characterized by different turbulence 
conditions. The optical turbulence has been calculated above these two sites for the 
same 15 nights studied for Dome C and a comparison between the three sites has been 
performed. 
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1 INTRODUCTION 

The Internal Antarctic Plateau represents a potential in- 
teresting location for astronomical applications. For almost 
a decade astronomers have shown more and more interest 
towards this region of the Earth thanks to its peculiar at- 
mospheric conditions. The extreme cold temperature, the 
dry atmosphere, the fact that the plateau is at more than 
2500 m above the sea level, that the turbulence seems to 
develop mainly in a thin surface layer of the order of 30-40 
m on the top of summits and that the seeing above this sur- 
face layer assumes values comparable to those obtained at 
mid-latitude sites get this region of the earth very appealing 
for astronomers. South Pole has been the first site equipped 
with an Observatory in the Internal Antarctic Plateau in 
which measurements of the optical turbulence have been 
done (Marks et al. 1996, Marks et al. 1999). Fifteen bal- 
loons have been launched in the winter period and it has 
been observed that the seeing above a surface layer of ~ 220 
m was very good (0.37 arcsec). Measurements of the optical 
turbulence at Dome C are more recent. After the first obser- 
vations done in 2004 with a MASS (Lawrence et al. 2004), 
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a series of studies done with different instrumentation have 
been published aiming to provide the assessment of the inte- 
grated seeing (Aristidi et al. 2005, Aristidi et al. 2009) and 
the vertical distribution of the optical turbulence (Trinquet 
et al. 2008). 

This paper deals with a different approach to the site 
assessment. In this context we are interested in investigat- 
ing the abilities of a mesoscale model (Meso-NH) in recon- 
structing correct optical turbulence features above different 
sites of the Internal Antarctic Plateau and its abilities in 
discriminating the optical turbulence properties of different 
sites. Meso-NH (Lafore et al., 1998) is a non- hydrostatic 
mesoscale research model developed jointly by the Centre 
National des Recherches Meteorologiques (CNRM) and the 
Laboratoire dAerologie de Toulouse, France. The Astro- 
Meso-NH package (Masciadri et al. 1999a) was first proven 
to be able to reconstruct realistic profiles above astro- 
nomical sites by Masciadri et al. (1999b) and Masciadri et 
al. (2001) and statistically validated later on (Masciadri & 
Jabouille, 2001, Masciadri et al. 2004, Masciadri & Egner 
2006). In the Astro-Meso-NH package all the main inte- 
grated astroclimatic parameters such as the isoplanatic an- 
gle, the wavefront coherence time, the scintillation rate, the 
spatial coherence outer scale are coded in the model (Mas- 
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Table 1. Results obtained in Lascaux et al. (2010a) that proved the Meso-NH model reliability above Dome C in reconstructing the 
optical turbulence spatial distribution. Three parameters are estimated: the mean surface turbulent layer (h s j), the seeing in the free 
atmosphere (efa) obtained integrating the from h s i up to the end of the atmosphere, the total seeing (stot) obtained integrating 
the from the ground up the top of the atmosphere. Beside each parameter is reported the associated standard deviation (cr) and the 
statistical error (a/VN). 







a 


a/VN 




a 


a/VN 


stot 


a 


a/VN 




(m) 






(arcsec) 






(arcsec) 






Observations 


35.3 


19.9 


5.1 


0.30 


0.70 


0.20 


1.60 


0.70 


0.20 


Model 


44.2 


24.6 


6.6 


0.30 


0.67 


0.17 


1.70 


0.77 


0.21 



Table 2. Geographic coordinates of Dome A, Dome C and South Pole. The altitude is in meter. 



SITE 


LATITUDE 


LONGITUDE 


MESO-NH 


MEASURED 








ALTITUDE (m) 


ALTITUDE (m) 


Dome A* 


80°22'00"S 


077°21'11"E 


4089 


4093 


Dome C** 


75°06'04"S 


123°20'48"E 


3230 


3233 


South Pole 


90°00'00"S 


000°00'00"E 


2746 


2835 



* GPS measurement by Dr. X. Cui (private communication). 

** GPS measurement by Prof. J. Storey (private communication). 



ciadri et al. 1999a). The model is also coded to calculate 
the astroclimatic parameters in finite vertical slabs (h m i n , 
hmax) in the troposphere (Masciadri et al. 1999a, Masciadri 
& Garfias 2001, Lascaux et al. 2010b). It can be therefore a 
useful tool for adaptive optics applications in classical as well 
as GLAO and/or MCAO configurations because we can pro- 
duce OT vertical distribution in whatever vertical slab we 
wish and with the suited vertical resolution. More recently, 
Meso-NH has been statistically validated above Dome C by 
Lascaux et al. (2009, 2010a). The most important results 
obtained in these two last papers are summarized in Table 
1. Briefly, the observations at Dome C, for a set of 15 win- 
ter nights (all the available nights for which is known the 
optical turbulence vertical distribution), gave a mean sur- 
face layer thickness h s ^ bs = 35.3 =b 5.1 m. The simulated 
surface layer thickness obtained with the Meso-NH model 
{hsi^mnh — 44.2 ± 6.6 m) is well correlated to measurements. 
The statistical error is of the order of 5-6 m but the stan- 
dard deviation (cr) is of the order of 20-25 m. This indicates 
that the statistic fluctuation of this parameter is intrinsically 
quite important. The median simulated free- atmosphere see- 
ing (s m nh,FA = 0.30 =b 0.17 arcsec) as well as the median 
total seeing (s m nh,TOT = 1.70 ± 0.21 arcsec) are well corre- 
lated to observations, respectively e bs,FA = 0.3 ± 0.2 arcsec 
and s bs,TOT = 1.6 ± 0.2 arcsec. 

In the context of this paper we consider that the Meso- 
NH model is calibrated as shown in Lascaux et al. (2010a) 
i.e. it produces optical turbulence features in agreement with 
observations. We therefore apply the Meso-NH model with 
the same configuration to other two sites of the plateau: 
South Pole and Dome A (Table 2). 

Why these sites ? Dome A is an almost uncontaminated 
site of the plateau. It is the highest summit of the plateau 
and, for this reason, it is expected to be among the best 
astronomical sites for astronomical applications. The high 
altitude reduces the whole atmospheric path for light coming 
from space and above the summit the katabatic wind speed 
is reduced to minima values. Dome A has been proved to 



Table 3. Meso-NH model configuration. In the second column 
the horizontal resolution AX, in the third column the number 
of grid points and in the fourth column the horizontal surface 
covered by the model domain. 



Domain 


AX 


Grid Points 


Surface 




(km) 




(kmxkm) 


Domain 1 


25 


120x120 


3000x3000 


Domain 2 


5 


80x80 


400x400 


Domain 3 


1 


80x80 


80x80 



have the strongest thermal stability (Hagelin et al. 2008) 
in proximity of the ground due to the coldest temperature. 
Dome A is a Chinese base. In the last few years the Chinese 
astronomers gave a great impulse to the site characterization 
showing a great interest for building astronomical facilities 
in this site. Optical turbulence measurements during the 
winter time are not yet available but site testing programs 
are on-going (Ashley et al. 2010). South Pole is interesting 
in our study because measurements of optical turbulence 
are available and, at the same time, the site is not located 
on a summit but on a gently slope. From the preliminary 
measurements done in the past we expect a surface turbulent 
layer that is thicker than the surface layer developed above 
the other two sites (Dome C and Dome A) due to the ground 
slope and the consequent katabatic winds in proximity of 
the surface. The three sites form therefore a perfect sample 
for a benchmark test on the model behavior and the model 
abilities. 

In Section 2 the numerical set-up of the model is pre- 
sented. In Section 3 results of the complete analysis of the 
three major parameters that characterize the optical tur- 
bulence features: surface layer thickness, seeing in the free 
atmosphere i.e. calculated above the surface layer and to- 
tal seeing are reported. Two different criteria to define the 
surface layer are used with consequent double treatment. Fi- 
nally, in Section 4 the results of this study are summarized. 
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2 NUMERICAL SET-UP 

Meso-NH (Lafore et al. 1998) can simulate the temporal 
evolution of the three-dimensional atmospheric flow over any 
part the globe. The prognostic variables forecasted by this 
model are the three cartesian components of the wind u, 
v, w, the dry potential temperature G, the pressure P, the 
turbulent kinetic energy TKE. 

The system of equation is based upon an anelastic 
formulation allowing for an effective filtering of acoustic 
waves. A Gal-Chen and Sommerville (1975) coordinate on 
the vertical and a C-grid in the formulation of Arakawa 
and Messinger (1976) for the spatial digitalization is used. 
The temporal scheme is an explicit three-time-level leap- 
frog scheme with a time filter (Asselin 1972). The turbulent 
scheme is a one-dimensional 1.5 closure scheme (Cuxart et 
al. 2000) with the Bougeault and LacarrereBougeault, P. 
& Lacarrere (1989) mixing length. The surface exchanges 
are computed in an externalized surface scheme (SUR- 
FEX) including different physical packages, among which 
ISBA (Noilhan et al. 1989) for vegetation. Masciadri et al. 
(1999a,b) implemented the optical turbulence package to be 
able to forecast also the optical turbulence {C% 3D maps) 
and all the astroclimatic parameters deduced from the C%. 
We will refer to the Astro-Meso-NH code' to indicate this 
package. The integrated astroclimatic parameters are calcu- 
lated integrating the C% with respect to the zenith in the 
Astro-Meso-NH code. We list here the main characteristics 
of the numerical configuration used in this study: 

• The interactive grid-nesting technique (Stein et al. 
2000) is used, with three imbricated domains of increased 
horizontal mesh-sizes (AX=25 km, 5 km and 1 km, Table 3). 
Such a method is used to permit us to achieve the best reso- 
lution on a small surface but keeping the volumetric domain 
in which the simulation is done in thermodynamic equilib- 
rium with the atmospherical circulation that evolves at large 
spatial scale on larger domains. We shown (Lascaux et al. 
2009) that the simulations results are sensitive to the chosen 
horizontal resolution. To achieve a good correlation between 
model outputs and observations, a grid-nesting configura- 
tion with a high horizontal resolution (at least AX=1 km) 
is mandatory. 

• The vertical grid is the same for all the domains re- 
ported in Table 3. The first vertical grid point is at 2 m 
above ground level (a.g.l.). A logarithmic stretched grid up 
to 3500 m a.g.l. (with 12 points in the first hundred of me- 
ters) is employed. Above 3500 m a.g.l., the vertical resolution 
is constant (AH ~ 600 m). The maximum altitude achieved 
is around 20 km a.g.l.. The first point at only 2 m above the 
ground (and with 12 points in the first hundred of meters) 
is necessary to forecast the typical very thin surface layer 
observed in the Ant art ic Plateau. 

• All simulations are initialized and forced every 6 hours 
at synoptic times (00:00, 06:00, 12:00, 18:00) UTC by anal- 
yses from the European Center for Medium-range Weather 
Forecasts (ECMWF) 1 . The simulations run for 18 hours. 
Note that the time at which the simulation starts (UTC) 
differs for Dome A, Dome C and South Pole. This is done so 
to be able to compare optical turbulence profiles simulated 

1 ECMWF: http : // www.ecmwf.int / 



in the same temporal interval with respect to the local time 
(LT). For each night, a mean vertical profile of C% is com- 
puted between the time interval (20:00 - 00:00) LT as done 
in Lascaux et al. (2009, 2010a). This range is centered on 
the time at which the balloons were typically launched at 
Dome C. In this way we obtain the most representative sim- 
ulated Cfj profile for each night 2 . In Table 4 are reported, 
for each site, the time at which the simulation starts and 
the duration AT of the simulation with respect to the local 
time. 

• An optimized version of the externalized surface scheme 
ISBA (Interaction Soil Biosphere Atmosphere) for antarctic 
conditions is employed (Le Moigne et al. 2009, 2010). Such 
a scheme has been used in Lascaux et al. (2010a) and it 
contributed to provide a realistic reconstruction of the opti- 
cal turbulence near the surface (optical turbulence strength 
and turbulence layer thickness) . It is indeed obvious that the 
most critical part of an atmospherical model for this kind 
of simulations is the scheme that controls the air/ground 
turbulent fluxes budget. Our ability in well reconstructing 
the surface temperature T s is related to the ability in recon- 
structing the sensible heat flux H that is responsible of the 
buoyancy-driven turbulence in the surface layer. 

• The Astro-Meso-NH package (Masciadri et al. 1999a) 
implemented in the most recent version of Meso-NH has 
been used to calculate the optical turbulence and derived 
astroclimatic parameters. 

As shown in Lascaux et al. (2010a), the best choice for 
the description of the orography is the RAMP (Radarsat 
Antarctic Mapping Project) Digital Elevation Model (DEM) 
presented in Liu et al. (2001), instead of the GTOPO30 
DEM from the U.S. Geological Survey used in Lascaux et 
al. (2009). For this study, therefore, the RAMP Digital El- 
evation Model has been used. The orography of each area 
of interest in this study (Dome C, Dome A, South Pole) is 
displayed on Fig. 1. All the grid-nested domains, from low 
horizontal resolution (larger mesh-size) to high horizontal 
resolution (smaller mesh-size) are reported. As can be seen 
in Fig. 1 (c,f,i) the orography around Dome C and Dome 
A is more detailed than the orography in proximity of the 
South Pole. This is due to the fact that the procedure to 
obtain a DEM integrates data from many different sources 
(satellite radar altimetry, airborne surveys, GPS surveys, 
station-based radar sounding...). However the resolution of 
some areas (typically those that can hardly receive infor- 
mation from the satellites) remain poorer than others. The 
region included in the inner circular polar region (and there- 
fore South Pole) fits with this condition and this is the rea- 
son why the orography is somehow less detailed than the 
rest of the Internal Antarctic Plateau. Nevertheless, this is 

2 In the prediction of a parameterized parameter (such as the op- 
tical turbulence) there is not a 1-1 correlation with the real time. 
This means, with an explicative example, that is somehow mean- 
ingless to predict the turbulence at a precise time t=t* as we do 
for a parameter that we resolve explicitly such as the temperature 
or the wind speed. This is the reason why, in order to obtain the 
most representative profile to be compared to measurements, 
we calculate the mean of the in a temporal interval AT. Such 
a procedure has been used in many previous papers (Masciadri 
and Jabouille 2001; Masciadri et al. 2004; Masciadri and Egner 
2006) 
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Table 4. Simulation starting time and time interval chosen for computations for the 3 different sets of simulations (Dome A, Dome 
C and South Pole). 





Dome A 


Dome C 


South Pole 


Starting time 


06:00 UTC / 11:00 LT 


00:00 UTC / 08:00 LT 


12:00 UTC / 12:00 LT 


Time interval 
for C^j computations 
(20:00 - 00:00 LT) 


15:00 - 19:00 UTC 


12:00 - 16:00 UTC 


20:00 - 00:00 UTC 



a region with no peaks or mountains and with just a regular 
and gently slope. We can therefore reasonably expect that 
the poorer accuracy in the orography has little or minor in- 
fluence on the results of the numerical simulations done with 
a mesoscale model such as Meso-NH. 

The same set of 15 winter nights used by Lascaux et al. 
(2009, 2010a) to validate the model above Dome C is inves- 
tigated in this study for the three antarctic sites Dome C, 
Dome A and South Pole. 



3 OPTICAL TURBULENCE ABOVE DOME C, 
DOME A AND SOUTH POLE 

In this section we investigate and compare the values ob- 
tained above the three sites (Dome C, Dome A and South 
Pole) of three parameters that characterize the optical tur- 
bulence features above the antarctic plateau: 

• surface layer thickness; 

• free atmosphere seeing from the surface layer thickness 
(h s i) up to the top of the atmosphere; 

• total seeing from the ground up to the top of the atmo- 
sphere. We note that this corresponds to ~10 km because 
the balloons explode at this altitude due to the high pressure 
and the strong wind speed. 

In a numerical mesoscale model the great challenge and diffi- 
culty is related to the parameterization of the optical turbu- 
lence. The critical issue is related to the ability of the model 
in reconstructing the vertical distribution of the optical tur- 
bulence (i.e. the C%). This is the reason why we selected and 
studied these three fundamental parameters. The other in- 
tegrated astroclimatic parameters are obtained calculating 
the integral of the C% and wind speed vertical profiles along 
the troposphere. A forth-coming paper will be dedicated to 
the analysis of the integrated astroclimatic parameters. 



3.1 Optical turbulence surface layer thickness 

To compute the surface layer thickness for each night, the 
same method employed in Trinquet et al. (2008) and Las- 
caux et al. (2010a) is first used. The thickness h s i is denned 
as the vertical slab containing 90 per cent of the optical 
turbulence developed inside the first kilometer above the 
ground: 



lkrr 
8m 



C 2 N (h)dh 



< 0.90 



(1) 



Table 5. Mean surface layer thicknesses h s i computed for the 3 
sites, for 15 different winter nights using the criterion in Eq. 1. 
Units in meter (m). The mean values are also reported with the 
associated statistical error cr/y/N. 



Date 


Dome A 


Dome C 


South Pole 


04/07/05 


65.0 


30.4 


117.6 


07/07/05 


529.4 


35.4 


262.9 


11/07/05 


28.6 


80.0 


131.9 


18/07/05 


27.7 


49.7 


224.0 


21/07/05 


17.6 


66.7 


136.3 


25/07/05 


15.7 


27.4 


298.6 


01/08/05 


25.5 


22.6 


185.2 


08/08/05 


53.1 


34.2 


104.4 


12/08/05 


19.4 


16.7 


59.0 


29/08/05 


17.4 


91.4 


251.3 


02/09/05 


16.4 


70.9 


164.9 


05/09/05 


125.2 


338.4 


128.0 


07/09/05 


59.8 


52.5 


103.6 


16/09/05 


38.8 


19.4 


158.7 


21/09/05 


20.1 


21.0 


148.0 


Mean 


37.9* 


44.2* 


165.0 


cr 


30.2* 


24.6* 


67.3 




8.1* 


6.6* 


17.4 



where is the refractive index structure parameter. We 
remind here that the selection of this criterium (that we call 



*These values are computed without taking into account the 
night of the 05/09/05 for Dome C and the night of 07/07/05 
for Dome A (see text for further explanations). 



criterium A) is motivated by the fact that we intend to com- 
pare our calculations with measurements done by Trinquet 
et al. (2008). This criterium has been selected by Trinquet 
et al. (2008) because the typical optical turbulence features 
above the internal antarctic plateau is characterized by a 
major bump at the surface and a consistent decreasing of 
the optical turbulence strength in the first tens of meters. 
The selection of the percentage is obviously absolutely ar- 
bitrary and, in this context, is mainly useful to check the 
correlation with measurements and to compare predictions 
on different sites (in relative terms therefore). The choice of 
the inferior limit of the integral (8 m) is motivated by the 
fact that Trinquet et al. (2008) intended to compare results 
obtained with balloons with those provided by the DIMM 
placed at 8 m from the ground. 

Table 5 reports the computed values of the surface layer 
thickness for each night at the three sites, as well as the 
mean, the standard deviation (a) and the statistical error 
(cr/y/N) for the 15 nights. For each night, the surface layer 
thickness is computed from a computed C% profile aver- 
aged between 20 LT and 00 LT (see Table 4 for hours in 
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Figure 1. Orography of three different regions of the internal Antarctic Plateau as seen by the Meson-Nh model (polar stereographic 
projection, grid-nesting configuration), (a), (b) and (c) show the three imbricated domains for the Dome C simulations, with horizontal 
resolution of 25 km, 5 km and 1 km, respectively, (d), (e) and (f) show the three imbricated domains for the Dome A simulations, with 
horizontal resolution of 25 km, 5 km and 1 km, respectively, (g), (h) and (i) show the three imbricated domains for the South Pole 
simulations, with horizontal resolution of 25 km, 5 km and 1 km, respectively. The dot labeled 'C indicates the Concordia Station. The 
dot labeled 'A' indicates Dome A. SP stands for South Pole. The altitude is expressed in meter (m). 



UT) as done in Lascaux et al. (2010a). The calculated mean 
surface layer thickness above South Pole is h s j = 165 m ± 
17.4 m, at Dome C h s j=44.2 m zb 6.6 m and at Dome A 
h s z=37.9 m zb 8.1 m. In this paper we are not forced any- 
more to use the same inferior limit of the integral in Eq. 1 
(8 m) than Trinquet et al. (2008), and we can compute the 
surface layer thickness starting the integral at the ground. 
Under this assumption the calculated mean surface layer 
thickness at South Pole is h s i =158.7 m zb 16.2 m, at Dome 
C h^=45.0 m ± 7.1 m and at Dome A h sZ =34.9 m zb 7.9 m. 
We conclude that at South Pole, h s i is more than three time 
larger than at Dome C or Dome A in both cases. This differ- 
ence is well correlated with previous observations done above 



South Pole. More precisely, observations related to 15 bal- 
loons launched during the period (20/6/1995 - 18/8/1995) 
indicated h s j=220 m (Marks et al. 1999). Measurements in 
that paper are done in winter but in a different year and for 
different nights. It is not surprising therefore that the match- 
ing between calculations and measurements is not perfect. 
Unfortunately the precise dates of nights studied in the pa- 
per from Marks et al. (1999) are not known. It is therefore 
not possible to provide a more careful estimate. It is how- 
ever remarkable that the h s i above South Pole is substan- 
tially larger than the h s i above Dome C and Dome A. Also 
we note that the typical thickness calculated above South 
Pole with a statistical sample of three months by Swain and 
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Gallee (2006) was h s ^=102 m. The authors used however a 
different definition of turbulent layer thickness. More pre- 
cisely, they defined h s i as the elevation (starting from the 
lowest model level) at which the turbulent kinetic energy 
contains 1 per cent of the turbulent kinetic energy of the 
lowest model layer. A comparison of this result with our cal- 
culations and with measurements is therefore meaningless. 
The same conclusion is valid for the estimates of h s i given at 
Dome C as already explained in Lascaux et al. (2009, 2010a). 
In conclusion, looking at Table 5, individuals values for each 
nights show a h s ^sp almost always higher than 100 m, with 
a maximum close to 300 m (2005 July 25), whereas Ji s i,dc 
and h s i,DA are always below 100 m. Dome C and Dome A 
have a comparable surface layer thickness. For this sample 
of 15 nights, h s i,DA is 6.3 m smaller than h s i,DC- We note 
also that the number of nights for which h s i is very small 
(inferior at 30 m) is more important at Dome A (nine in- 
stead of six at Dome C). This difference is however not really 
statistically reliable considering the number of the nights in 
the sample. For a more detailed discrimination between the 
h s i value at Dome C and Dome A we need a larger statistic. 
This analysis is planned for a forthcoming paper. 

Looking at the results obtained night by night we can 
note some specific features observed in specific cases. Two 
nights (September 5 at Dome C (h s i = 338.4 m) and July 
7 at Dome A (h s i = 529 .4 m)) present similar characteris- 
tics: the surface layer thickness h s i is well larger than the 
observed one. In these two cases, however, as already ex- 
plained in Lascaux et al. (2009, 2010a) for the case of Dome 
C, the large value of h s i does not mean that a thicker and 
more developed turbulence is present near the ground but it 
simply means that, in the first kilometer from the ground, 90 
per cent of the turbulence develops in the (0, h s i) range. On 
September 5, at Dome C, the model reconstructs the total 
seeing on the whole 20 km much weaker than what has been 
observed and more uniformly distributed and, consequently, 
the criterium (Eq.l) provides us a much larger value of h s i. 
In both cases (on September 5 at Dome C and on July 7 at 
Dome A), when we look at the vertical distribution of the 
C]^ calculated by the model, we observe that the turbulence 
is concentrated well below 20 m in a very thin surface layer 
with a very weak total seeing (see next section). The case 
of 5 September at Dome A, is however a case in which the 
model reconstructed a surface turbulent layer thicker than 
what has been observed. 

It is known that mesoscale model provide a temporal 
variability of the turbulence in the high part of the atmo- 
sphere that is smoother than what observed with vertical 
profiler. This is due to the fact that a mesoscale model is 
more active in the low part of the atmosphere where the 
orographic effects are mainly present. We recently obtained 
(Lascaux et al. 2009) very encouraging results showing that 
the C% in the free atmosphere has a temporal variability 
even on a small dynamic range (-18, -16.5 in logarithmic 
scale). This is a signature of the improvement of the model 
activity in the high part of the atmosphere. At present, how- 
ever, it presents a hazard to quantifying the typical time- 
scale for temporal variability of all the parameters related to 
the optical turbulence reconstructed by a mesoscale model. 
Nevertheless we can describe the temporal variability of the 
morphology of these parameters such as, for example the 
thickness of the surface layer. 



In Appendix A, we report the temporal evolution of 
the calculated C% for all the nights above the three sites. 
Looking at these pictures we can give a description of the 
morphology of the temporal variation of the surface layer. 
Above Dome C and Dome A, the thickness of the surface 
layer remains mostly stable during the night, even though 
we have night to night variations as shown in Table 5. This 
fits with preliminary results shown in Ashley et al. (2010) 
above Dome A. Above South Pole, the thickness varies in a 
much more important way during the night with oscillations 
that can reach 50 to 100 m. The larger variability of the 
typical turbulent surface layer thickness is confirmed also 
by the larger value of a observed above South Pole (Table 
5,6,7). 

In order to compare our calculations and results with 
those obtained by Swain and Gallee (2006) we applied also 
a different criterium (criterium B) based on the analysis of 
the vertical profile of turbulent kinetic energy (TKE) in- 
stead of the vertical profile of Cfj. The TKE is certainly 
an ingredient from which the optical turbulence depends on 
and it represents the dynamic turbulent energy. However, 
it is known (Masciadri & Jabouille 2001) that the C% de- 
pends also on the gradient of the potential temperature and 
moreover, the selection of the value of percentage of the tur- 
bulent kinetic energy (1%, 10%, other...) used as a threshold 
is absolutely arbitrary. This method is therefore not useful 
to quantify the absolute value of h s i to be compared to mea- 
surements provided by Trinquet et al. (2008) and Marks et 
al. (1999). It can possibly be useful for relative comparisons 
between different sites or to compare our calculations with 
calculations provided by Swain Sz Gallee (2006). 

Using this method (Table 6), the surface layer height 
is determined as the elevation at which the TKE is X% of 
the lowest elevation value. We calculated the h s i for X = 1 
(Table 6) and X = 10 (Table 7). X=l is the case treated by 
Swain and Gallee (2006). For each simulation, we first com- 
pute the average of the TKE profile for the night between 
20 LT and 00 LT. While the average of the C 2 N profile is 
calculated with a 2 minutes rate sample, the average of the 
TKE is calculated with 5 profiles, available at each hour (20, 
21, 22, 23, 00) LT. This gives us an averaged vertical profile 
of TKE characteristic of the considered night. 

The computation of the surface layer thickness is then 
performed using this averaged TKE profile. It has been ob- 
served that, when the night presents only low dynamic tur- 
bulence (with a very low averaged TKE at the lowest eleva- 
tion level), it is very hard to retrieve a surface layer height 
using this criterium. This means that the turbulence is so 
weak that we are at the limit of necessary turbulent kinetic 
energy to resolve the turbulence itself. For these nights (in- 
dicated with an asterisk in Table 6) it could happen that 
we calculated the average on a number of estimates smaller 
than 5 (as for all the other cases). The results are reported 
in Table 6. 

Table 6 shows that results obtained with the criterium 
of the TKE are similar to those obtained with the criterium 
described in Eq.l. Table 7 provides smaller values of h s i 
above all the three sites. We treat the case (X = 10) to 
show that, tuning the value of the percentage, it is possible 
to find different values of h s i. This means that h s i estimates 
are useful only if they are compared to measurements using 
the same criteria. 



Optical turbulence at Dome C, Dome A and South Pole 7 



Table 6. Mean surface layer thicknesses h s i computed for the 3 
sites, for the same set of nights shown in Table 5, but computed 
with a different criterion. The surface layer height is determined 
as the elevation at which the averaged TKE between 20 LT and 
00 LT for each nigh is 1% of the averaged lowest elevation value. 
Units in meter (m). The mean values are also reported with the 
associated statistical error a/y/N. 



Date 


Dome A 


Dome C 


South Pole 


04/07/05 


78* 


32 


112 


Ul / U 1 1 uo 


D 


oz 


1 1 9* 


11/07/05 


40 


76 


174 


18/07/05 


32 


48 


242 


21/07/05 


22 


56 


144 


25/07/05 


22 


12* 


148 


01/08/05 


32 


22 


186 


08/08/05 


56 


32 


112 


12/08/05 


22 


60 


58 


29/08/05 


22 


82 


250 


02/09/05 


20 


60 


188 


05/09/05 


136 


30* 


146 


07/09/05 


72 


74 


250 


16/09/05 


56 


22 


192 


21/09/05 


26 


22 


170 


Mean 


42.8 


44 


165.6 


a 


33.0 


22.6 


55.5 


a/y/N 


8.5 


5.8 


14.3 



*These values are computed using a number 
of profiles minor than five. 



Table 7. Mean surface layer thicknesses h s i computed for the 3 
sites, for the same set of nights shown in Table 5, but computed 
with a different criterion. The surface layer height is determined 
as the elevation at which the averaged TKE between 20 LT and 
00 LT for each nigh is 10% of the averaged lowest elevation value. 
Units in meter (m). The mean values are also reported with the 
associated statistical error a/y/N. 



Date 


Dome A 


Dome C 


South Pole 


04/07/05 


58 


22 


56 


07/07/05 


2 


24 


62 


11/07/05 


28 


52 


64 


18/07/05 


22 


32 


186 


21/07/05 


14 


40 


110 


25/07/05 


12 


6 


102 


01/08/05 


22 


16 


126 


08/08/05 


40 


24 


102 


12/08/05 


16 


8 


40 


29/08/05 


14 


68 


192 


02/09/05 


14 


46 


112 


05/09/05 


106 


30 


88 


07/09/05 


50 


52 


68 


16/09/05 


40 


12 


150 


21/09/05 


18 


14 


116 


Mean 


30.4 


27 


104.9 


a 


26.1 


15.6 


45.2 


a/y/N 


6.7 


4 


11.7 



To conclude, both criteria (A and B with X=l) give sim- 
ilar mean h s i values for all the 3 sites for this limited set of 
nights. Evaluating the surface layer thickness over a more ex- 
tended set of nights should be the next step. It would permit 
us to compute more reliable and robust statistical estimates 
for h s i over the 3 antarctic sites and possibly to identify dis- 
crimination between the h s i at Dome A and Dome C. This 
also means that our estimate of h s j = 165 m above South Pole 
is better correlated to measurements (h s z=220 m) than the 
estimate (h s ^ = 102 m) obtained by Swain and Gallee (2006) 
at the same site. 



3.2 Seeing in the free atmosphere and seeing in 
the whole atmosphere 

The seeing in the free atmosphere and in the whole atmo- 
sphere for A=0.5xl0 _6 m is: 

/ rhtop \ 3 / 5 

efa = 5.41 • A" 1/5 • C 2 N (h)-dhj (2) 

£ TOT = 5.41 • A" 1/5 • ( J C 2 N (h)-dh) (3) 

with h top ~ 13 km from the sea level i.e. where the balloons 
explode and we have no more their signal. Table 8 shows 
the simulated total seeing (eror) and free- atmosphere see- 
ing (sfa) for each night and each sites (Dome C, Dome A 
and South Pole). We define the free atmosphere as the por- 
tion of the atmosphere extended from the mean h s i reported 
in Table 5 up to htop- The median values of the seeing as 
well as the standard deviation (a) and the statistical er- 
ror (a/y/N) are reported. As expected the total seeing is 
stronger at Dome A (stot,da = 2.37 =b 0.27 arcsec) than 
at Dome C (etot,dc = 1.70 ± 0.21 arcsec) or South Pole 
(stot,sp = 1.82 ± 0.23). The total seeing is very well corre- 
lated with measurements at Dome C (Lascaux et al. 2010a - 
£roT,obs = 1.6 arcsec) and at South Pole (Marks et al. 1999 
- 6TOT,obs = 1-86 arsec) getting the estimate at Dome A 
highly reliable. The minimum median free- atmosphere see- 
ing is found at Dome A (efa,da = 0.23 =b 0.28 arcsec). 
The medain free- atmosphere seeing at Dome C is 6fa,dc 
= 0.30 ± 0.17 arcsec and at South Pole, s F a,sp = 0.36 ± 
0.11 arcsec. The seeing in the free atmosphere is very well 
correlated with measurements at Dome C (Lascaux et al. 
2010a - SFA,obs = 0.30 arcsec) and at South Pole (Marks 
et al. 1999 - 6FA,obs = 0.37 arsec) getting again very reli- 
able the method (Meso-NH model) as well as the estimates 
at Dome A. What is remarkable is that, even if h s i,DA < 
h s i, dc < h s i,sp. Dome A is the site with the lowest free- 
atmosphere seeing efa- This means that at Dome A as well 
as at Dome C the turbulence is concentrated inside the first 
tens of meters from the ground. Moreover, the turbulence 
in the surface layer is stronger at Dome A than at Dome C. 
This can be explained with the stronger thermal stability 
of Dome A near the ground. Our results match, therefore, 
with predictions we did in Hagelin et al. (2008) studying 
only features of the meteorological parameters. 

At South Pole, however, the C% vertical distribution 
decreases in a less abrupt way because the thermal stability 
near the ground is less important. The vertical distri- 
bution is spread over hundreds of meters from the ground, 
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Table 8. Total seeing £TOT=£[8m,h t o P ] an( ^ see i n g hi the free 
atmosphere £FA=£[h s i,h top ] calculated for the 15 nights and av- 
eraged in the temporal range 20-00 LT. See the text for the defi- 
nition of h s i and ht op - 





DOME A 


DOME C 


SOUTH POLE 


Date 


efaI^tot 


efaI^tot 


efa/ £tot 




(h sZ =37.9m) 


(h sZ =44.2m) 


(h sZ =165m) 


04/07/05 


2.55 / 3.37 


0.22 / 2.28 


0.40 / 1.67 


D7 /D7 /n^i 


n 9n / n 94 

U.ZU / (J.Zi^-t 


n or / i Qi 


U.O-L / U./U 


11/07/05 


0.23 / 2.78 


1.61 / 1.81 


0.47 / 1.96 


18/07/05 


0.21 / 2.73 


0.80 / 1.94 


1.46 / 2.28 


21/07/05 


0.21 / 1.95 


0.86 / 1.27 


0.31 / 1.71 


25/07/05 


0.22 / 1.55 


0.25 / 0.85 


0.32 / 0.76 


01/08/05 


0.22 / 1.78 


0.22 / 2.27 


0.52 / 1.78 


08/08/05 


1.45 / 2.42 


0.35 / 1.70 


0.28 / 1.69 


12/08/05 


0.23 / 2.37 


0.23 / 0.99 


0.29 / 1.82 


29/08/05 


0.23 / 1.83 


2.29 / 2.47 


1.55 / 2.11 


02/09/05 


0.22 / 1.76 


1.16 / 1.54 


0.81 / 3.56 


05/09/05 


3.21 / 3.36 


0.30 / 0.52 


0.31 / 2.98 


07/09/05 


2.43 / 3.49 


1.69 / 3.73 


0.31 / 1.41 


16/09/05 


1.11 / 4.60 


0.21 / 1.57 


0.99 / 3.96 


21/09/05 


0.20 / 2.30 


0.26 / 1.63 


0.36 / 2.32 


Median 


0.23 / 2.37 


0.30 / 1.70 


0.36 / 1.82 


a 


1.08 / 1.03 


0.67 / 0.77 


0.43 / 0.90 


a/VN 


0.28 / 0.27 


0.17 / 0.21 


0.11 / 0.23 



instead of tens of meters like for Dome A or Dome C. As a 
consequence the total seeing is also weaker than above Dome 
C and Dome A. 

Such a behavior is evidenced in Figure 2, which displays 
the median vertical profiles over the 3 sites. 

Looking at Table 8, we note that the values of a for 
the total seeing above the three sites is mostly comparable 
with no significant differences even if Dome C seems a little 
smaller (0.77) than Dome A (1.03) and South Pole (0.90). 
This indicates a comparable variability of the turbulence 
above the three sites. For what concerns the seeing in the 
free atmosphere, the value of a above Dome A is almost 
double (1.08) than above Dome C (0.67) and South Pole 
(0.43). 



4 CONCLUSION 

In this study the mesoscale model Meso-NH was used to 
perform forecasts of optical turbulence (evolutions of C% 
profiles) for 15 winter nights at three different antarctic 
sites: Dome A, Dome C and South Pole. The model has 
been used with the same configuration previously validated 
at Dome C (Lascaux et al. 2010a) and simulations of the 
same 15 nights have been performed above the three sites. 
The idea behind our approach is that once validated above 
Dome C, the model can be used above two other sites of 
the internal antarctic plateau to discriminate optical turbu- 
lence features typical of other sites. This should show the 
potentiality of the numerical tool in the context of the site 
selection and characterization in astronomy. South Pole has 
been chosen because in the past some measurements of the 
optical turbulence have been done and this can represent 



a useful constraint for the model itself. For Dome A there 
are not at present time measurements of the optical turbu- 
lence and this study provides therefore the first estimates 
ever done of the optical turbulence above this site. We test 
this approach above the antarctic plateau because this re- 
gion is particularly simple from the topographic point of 
view and certainly simpler than typical mid-latitude astro- 
nomical sites. No major mountain chains are present and the 
local surface circulations is mainly addressed by the energy 
budget air/ground transfer, the polar vortex circulation at 
synoptic scale and the katabatic winds generated by gravity 
effects on gently slopes due to the cold temperature of the 
iced surface. The main results we obtained are summarized 
in Table 9 and listed here: 

• We provide the first estimate of the optical turbulence 
extended on the whole 20 km above the Internal Antarctic 
Plateau. 

• The Meso-NH model achieves to reconstruct the three 
most important parameters used to characterize the optical 
turbulence: the turbulent surface layer thickness, the seeing 
in the free atmosphere and in the surface layer for the three 
selected sites: Dome C, Dome A and South Pole showing 
results in agreement with expectations. Measurements taken 
at Dome C and South Pole corresponds to balloons launched 
during 15 nights, in both cases. The statistic is not very large 
but reliable for a first significant result. The selected nights 
correspond to the 15 nights for which measurements of the 
Dome C are available. 

• Dome C and Dome A present a very thin surface layer 
size (h s i,DA = 37.9 =b 8.1 m and h s i,DC = 44.2 ± 6.6 m) 
while South Pole surface layer is much thicker (h s i,sp = 
165 =b 17.4 m). If we apply the criterium (A) described by 
Eq.(l) integrating from the ground instead than 8 m from 
the ground we find similar result within a couple of meters. 
All these estimates are well correlated with measurements. 
Surface layers calculated by the model at Dome C and Dome 
A have a comparable thickness considering the actual sam- 
ple. To better discriminate between the Dome A and Dome 
C surface layer thickness a richer statistic is necessary. An 
on-going study has started addressing this issue. 

• Dome A is the site with the strongest total seeing (2.37 
± 0.27 arcsec) with respect to Dome C (stot,dc = 1-70 =b 
0.21 arcsec) and South Pole (stot,sp = 1.82 ± 0.23 arcsec). 
This is explained by the stronger thermal stability near the 
ground with respect to the other two sites that cause large 
values of the optical turbulence in the thin surface layer. 

• All the three sites show a very weak seeing in the free at- 
mosphere i.e. above the correspondent mean h s i: 6fa,da = 
0.23 =b 0.28 arcsec at Dome A, efa,dc = 0.30 =b 0.17 arcsec 
at Dome C and £fa,sp = 0.36 ±0.11 arcsec at South Pole. 
Dome A show the weakest seeing in the free atmosphere. 

• The temporal variability of the thickness of the surface 
layer is more important at South Pole than above Dome A 
and Dome C that show a very stable trends in agreement 
with observations. The temporal variability of the seeing in 
the whole atmosphere does not show important differences 
above the three sites, while the variability of the seeing in 
the free atmosphere is almost double at Dome A than at 
Dome C and South Pole. 

• Both, the total seeing and the seeing in the free at- 
mosphere calculated by Meso-NH, are very well correlated 
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Figure 2. Median profiles simulated with the Meso-NH mesoscale model at Dome C (black), Dome A (blue) and South Pole (red). 
Left: from the ground up to 20 km. Middle: from the ground up to 1 km. Right: from the ground up to 200 m. Units are ra -2 / 3 . 



Table 9. Summary of the main results obtained in this study: surface layer h s ^, seeing in the free atmosphere (efa) an d total seeing 
(etot) at Dome C, Dome A and South Pole. The associated standard deviation (<x) and the statistical error (cr/y/N) are also reported. 



h s i a a/\/ r N £fa & cr/V^/V £tot & cr/V^/V 

(m) (arcsec) (arcsec) 



Observations - Dome C 


35.3 


19.9 


5.1 


0.30 


0.70 


0.20 


1.60 


0.70 


0.20 


Meso-NH - Dome C 


44.2 


24.6 


6.6 


0.30 


0.67 


0.17 


1.70 


0.77 


0.21 


Meso-NH - Dome A 


37.9 


30.2 


8.1 


0.23 


1.08 


0.28 


2.37 


1.03 


0.27 


Meso-NH - South Pole 


165.0 


67.3 


17.4 


0.36 


0.43 


0.11 


1.82 


0.90 


0.23 



with measurements at Dome C and South Pole getting the 
predictions done at Dome A highly reliable. 

• Dealing with the criteria used to define the surface layer 
thickness, we proved that, at least on the sample of 15 nights 
investigated, the criterium denned by Eq.l (criterium A) 
and the criterium using the vertical profile of the turbulent 
kinetic energy (TKE) taking h s i as the height at which the 
value of the TKE is less than 1% of the TKE at the low- 
est level near the ground (criterium B) provide very similar 
results. 

• The mean h s i we estimate at Dome C (h s j=44.2) is 
slightly thicker than what found by Swain and Gallee (2006) 
(h s z=27.7 m) with comparable discrepancy from measure- 
ments (h s i = 35.3 =b 5.1 m). The h s i we estimate at South 
Pole (h s i,sp = 165 =b 17.4 m) is thicker than what estimated 
by Swain and Gallee (2006) (h s i,sp = 102) but better cor- 
related to measurements (h s i,sp = 220 m) than what found 
by Swain and Gallee (2006). The h s i we estimate at Dome 
A (h s i, da — 37.9 =t 8.1 m) is somehow thicker than what 
estimated by Swain and Gallee (2006) (H s i,da = 18 m). It 
is however important to note that the standard deviation of 
h s i is of the order of h s i itself or even larger. The statistic 
error a/ \/~(N) is of the order of ~ 10 m. We think therefore 
that at present there are no major differences in our results 
with respect to Swain and Gallee (2006) with exeption of 
the fact that we proved that, with our model, the horizontal 
resolution of 1 km provides better results than a resolution 
of 100 km that is used by Swain and Gallee (2006). 

All these results deserve now a confirmation provided by an 
analysis done with a richer statistical sample. Also it would 
be interesting to refine this study when OT measurements 
above Dome A will be published. Besides, we can state that 
all major expectations concerning the typical features of the 



optical turbulence above South Pole, Dome C and Dome A 
have been confirmed by this study. The tendency shown by 
the model is obviously that in summer time, in proximity 
of the surface, due to the less stable regime, the turbulence 
thickness increases but the turbulence strength decreases. 
This is however, out the goals of this paper. 
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APPENDIX A: COMPUTED TEMPORAL 
EVOLUTIONS OF C% VERTICAL PROFILES 
FOR EACH NIGHTS AT DOME C, DOME A 
AND SOUTH POLE 

In this appendix we present all the individual figures of the 
18- hours temporal evolution of the C% for every night and at 
the three antarctic sites considered in this study (Figure Al: 
Dome A, Figure A2: Dome C and Figure A3: South Pole). 
The first couple of hours can be considered as spurious values 
because of the model adaptation to the ground. 
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Figure Al. Meso-NH temporal evolutions of vertical profiles at Dome A (log units) for the 15 forecasted nights, for all the 18 hours 
of the simulations, from the ground up to 400 m above ground level. Units are ra -2 / 3 . 
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Figure A2. Same than Figure Al but for Dome C. 
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Figure A3. Same than Figure Al but for South Pole. 



